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ROTATION ALLY SYMMETRIC BIHARMONIC MAPS BETWEEN MODELS 


S. MONTALDO, C. ONICIUC, AND A. RATTO 


Abstract. The main aim of this paper is to study existence and stability properties of rotationally 
symmetric proper biharmonic maps between two m-dimensional models (in the sense of Greene and 
Wu). We obtain a complete classification of rotationally symmetric, proper biharmonic conformal 
diffeomorphisms in the special case that m = 4 and the models have constant sectional curvature. 
Then, by introducing the Hamiltonian associated to this problem, we also obtain a complete de¬ 
scription of conformal proper biharmonic solutions in the case that the domain model is R 4 . In 
the second part of the paper we carry out a stability study with respect to equivariant variations 
(equivariant stability). In particular, we prove that: (i) the inverse of the stereographic projection 
from the open 4-dimensional Euclidean ball to the hyperbolic space is equivariant stable; (ii) the 
inverse of the stereographic projection from the closed 4-dimensional Euclidean ball to the sphere 
is equivariant stable with respect to variations which preserve the boundary data. 


1. Introduction 


Harmonic maps are critical points of the energy functional 

(1.1) E{ip) = \ [ \dip\ 2 dv g , 

1 Jm 

where ip : (M, g ) — > ( N , h ) is a smooth map between two Riemannian manifolds M and N. In 
analytical terms, the condition of harmonicity is equivalent to the fact that the map ip is a solution 
of the Euler-Lagrange equation associated to the energy functional (11.11) . i.e. 

(1.2) trace Vd(p = 0 . 

The left member of (11.21) is a vector field along the map tp, or, equivalently, a section of the pull¬ 
back bundle yW 1 ( TN ): it is called tension field and denoted r(^). In local charts, the tension field 
is given by the following expression: 

(1-3) T'ffoO = g ij (V(dv>))7. , 


where 

(1.4) 


(S7(dip))]j 


dV _ Mpfc , Ap7 ^ 

dx l dxi dx k a @ dx l dxi 


In (11.41) . M r and denote the Christoffel symbols of the Levi-Civita connections of ( M,g ) and 
(A, h) respectively. Also, note that Einstein’s convention of sum over repeated indices is adopted. 
We refer to mmmm for notation and background on harmonic maps. 
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A related topic of growing interest deals with the study of the so-called biharmonic maps: these 
maps, which provide a natural generalisation of harmonic maps, are the critical points of the 
bienergy functional (as suggested by Eells-Lemaire 0) 

(1.5) E 2 (p) = ^[ |t(<^)| 2 dv g . 

A Jm 

In p21J Jiang derived the first variation and the second variation formulas for the bienergy. In 
particular, he showed that the Euler-Lagrange equation associated to E 2 {p) is 

(1.6) r 2 (p) = — J {t{p)) = ~Xr(p) ~ traceR N (dy?, T(p))dp = 0 , 

where J is (formally) the Jacobi operator of ip, A is the rough Laplacian defined on sections of 
p- 1 ( TN) and 

(1.7) R N {X , Y) = V.vV Y ~ VyVx - V [x ,y\ 

is the curvature operator on ( N,h ). Equation (|1.6I) is a fourth order semi-linear elliptic system 
of differential equations. We also note that any harmonic map is an absolute minimum of the 
bienergy, and so it is trivially biharmonic; thus, a general working plan is to study the existence 
of proper biharmonic maps, i.e., biharmonic maps which are not harmonic. We refer to [T7] for 
existence results and general properties of biharmonic maps. 

Equivariant theory deals with special families of maps having enough symmetries to guarantee 
that harmonicity reduces to the study of a second order ordinary differential equation (we refer to 
m m noi nHusni [22] for background and examples). In [18], we developped a systematic approach to 
equivariant theory for biharmonic maps (see m, for recent developments). In this framework, the 
aim of this work is to study rotationally symmetric biharmonic maps between two m-dimensional 
models (in the sense of [II]). The study of rotationally symmetric biharmonic maps has been 
started in [3] and, in the particular case of maps between surfaces, in [21] , 

In this paper we obtain a complete classification of proper biharmonic conformal diffeomorphisms 
between two m-dimensional models in the special case that m = 4 and the models have constant 
sectional curvature (Proposition 13.11) . Then, by introducing the Hamiltonian associated to this 
problem, we obtain a complete classification of conformal proper biharmonic maps in the case that 
the domain model is M 4 and the target is any 4-dimensional model (Theorem 14. ip . 

In the second part of the paper we carry out a stability study with respect to equivariant vari¬ 
ations (equivariant stability). In particular, we prove that: (i) the inverse of the stereographic 
projection from the open 4-dimensional Euclidean ball to the hyperbolic space is equivariant sta¬ 
ble (Theorem 15.31) : (ii) the inverse of the stereographic projection from the closed 4-dimensional 
Euclidean ball to the sphere is equivariant stable with respect to variations which preserve the 
boundary data (Theorem 15.41) . 

Remark 1.1. We stated above that equations (11.21) and (11.61) are the Euler-Lagrange equations 
associated with the energy functional (jl.ll) and the bienergy functional (11.51) respectively. We point 
out that, in the case that the domain manifold M is not compact, the previous claim can be formally 
verified by considering these functionals over compact domains D C M, with smooth boundary 
dD , and smooth variations pt which coincide with p on M \int(U). In this context, a map p is a 
critical point if it is such over any compact domain D in M (see [2| for details). 

Acknowledgement. The authors wish to thank the referee for some useful comments that 
have improved the quality of the paper and, in particular, for pointing out the solutions given in 
Remark 13.61 
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2. Notation and preliminaries 

By definition, a point o of a Riemannian manifold M is said to be a pole if the exponential map 
exp 0 : T a M — > M is a diffeomorphism. Following mi, we say that an m-dimensional manifold 
(M m (o), g) with a pole o is a model if and only if every linear isometry of T 0 M can be realized 
as the differential at o of an isometry of M . The starting goal of Greene and Wu was to develop 
a good function theory on Cartan-Hadamard manifolds (CH-manifolds), i.e., complete simply- 
connected Riemannian manifolds of nonpositive sectional curvature. Models are more general than 
CH-manifolds (for instance, their curvature can be of variable sign) and turned out to be a suitable 
environment for Laplacian and Hessian comparison theorems. A significant geometric property of 
a model is the fact that we can describe it by means of geodesic polar coordinates centered at the 
pole o, as follows: 

(2.1) (M m (o), g ) = (S"*- 1 x [0, Too), f 2 (r)g s m -i + dr 2 ) , 
where 

(S m -\g S m- 1) 

denotes the (m — l)-dimensional Euclidean unit sphere, and the function /(r) is a smooth function 
which satisfies 

(2.2) /(0) = 0 , /'(0) = 1 and f(r) >0 if r > 0 . 

We also note that r measures the geodesic distance from the pole o. To shorten notation and 
emphasize the role of the function /, we shall write M™(o) to denote a model as in (12.11) . 

Remark 2.1. We observe that, if f(r) = r , then M™(o) = M m . We also point out that, if 
/(r) = (1/c) sinh(rc) (c > 0), then M™(o) represents H m (—c 2 ), i.e., the m-dimensional hyperbolic 
space of constant sectional curvature — c 2 . 

For future use (see |11|). we also recall that the radial curvature K(r) (r > 0) of a model Mj. n (o) 
is defined as the sectional curvature of any plane which contains d / dr . The radial curvature is 
related to the function f(r) by means of the following fundamental equation (the Jacobi equation): 

(2.3) /» T K(r) f (r) = 0 , /(0) = 0 , /'(0) = 1 . 

The aim of this paper is to study rotationally symmetric maps between two models, i.e. maps 
of the following type: 

Ta : Mf{p) -T M' m (o') 

(2.4) (0, r) i-t (6, a(r)) , 

where the function a(r) is smooth on [0, Too), positive on (0, Too) and, to ensure continuity, 
satisfies the boundary condition 

(2.5) a(0) = 0 . 

Remark 2.2. With a slight abuse of terminology, in our study of rotationally symmetric biharmonic 
maps we shall also consider the case that f{r) is defined on a finite interval [0, b], with f(b) = 0 and 
f'(b) = —1. In particular, we shall pay special attention to the case /(r) = (1/d) sin(rd), where 
d > 0 and 0 < r < (ir /d) : in this case, our manifold is the Euclidean m-sphere S m (d 2 ) of constant 
sectional curvature d?. 

Remark 2.3. More generally, one could apply the reduction technique of this paper to the following 
class of equivariant maps: 

: M} n (o ) -A M' h n (o') 

(6, r) i-t (W A (0), a(r)) , 
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( 2 . 6 ) 


where 1^(0) is a so-called eigenmap of eigenvalue A . That means that Ta : is a 

harmonic map with constant energy density equal to (A/2). Important examples of eigenmaps are: 
the identity map of (A = m — 1), the k-fold rotation e l8 e lkd of S 1 (A = k 2 ); and, also, 

the Hopf fibrations S ' 3 —>■ S 2 , S 7 —> S ' 4 and S ' 15 -» S 8 , with A equal to 8 , 16 and 32 respectively. 
However, since the primary aim of this paper is the study of conformal diffeomorphisms, we decided 
to focus on (12.41) . i.e., on the case that Ta( 0) is the identity map. 


By way of summary, the unknown function a(r) in (12.41) has to be determined in such a way 
that <p a be a biharmonic map. To this purpose, a calculation based on (11.41) shows that the tension 
field of a rotationally symmetric map as in (|2.4I) is given by 

d 


(2.7) 


T&a) = 


... . . .. f(r) h(a)h'(a ) 

a(r) + (m — 1) ! a(r) — (m — 1) 


da 


f(r ) / 2 (r) 

Therefore, in this case the reduced bienergy is given, up to an irrelevant constant factor, by the 
following expression: 


^ r+oo 

(2.8) E 2 {ipa) = 9 y o 


.. . . , . /'(r) • / \ / 1 \^'( a )^ , ( Q; ) 

a r) + (m - 1 77 -y a(r) - (m - 1 ) L \ 

Jv) f 2 yr) 


f m -\r)dr . 


We observe that the reduced bienergy functional (j2.8j) is of the form 


(2.9) 


E2(<Pa) = 


r*+oo 


L(r, a, a, a) dr , 


for a suitable Lagrangian function L. We showed in m that, in this symmetric context, the 
condition of biharmonicity is 

dL _ d_ f dL\ <P_ /d 
da dr \da) + dr 2 \ da y 

Now, a straightforward computation leads us to the explicit expression of (12.101) in our context (to 
simplify notation, we write / and h(a) instead of f(r) and h(a(r )) respectively): 


( 2 . 10 ) 


= 0 . 


f m 5 ^(m — l)h(a)(2ff"h'(a) — 2 (m — 3 )ffah"(a) + 2 (m — 4)f' 2 h'{a) 

— f 2 (h^ 3 \a)a 2 + 2 ah" (a)) + (m — 1 )h'(a) 3 ) + f((m — 3 ){m — 1 )f f' 2 a 

(2.11) — (m — 3 ){m — 1 )f' 3 d + (m — l)/(/(/^ 3 ^d + 2 f"a) — 2 ah'(a) 2 

— 3 a 2 h'(a)h" (a)) + ( m — l)f'(a((m — 4 )f f" — 2 (m — 3)h\a) 2 ) + 2f 2 a^ 3 ' > ) 
+ f 3 a ( 4 )) + (m — 1 ) 2 h(a) 2 h 1 (a)h"(a)^ = 0 . 


Remark 2.4. We point out that, to the purpose of comparison with the equation given, for m = 2, 
in (21l Corollary 2.3], equation (12.111) can be rewritten as follows: 


F" + (m - 1) 


ff'F' - h\a) 2 F 


j': 


p 


— (m — 1 ) 
h(a)h'(a) 


h{a)h" {a)F 

T 2 


= 0 


F = a + (m — 1)—d — (m — 1) 


In summary, a rotationally symmetric map ip a as in (12.41) is biharmonic if and only if a is a 
solution of ( 12 . 111 ) . 

In some cases it is useful to associate to a Lagrangian L as in (12.91) its corresponding Hamiltonian 
H, which is defined as follows (see, for instance, my- 


dL d dL 
da dr da 




( 2 . 12 ) 


H(r, a, a, a) = a 


4 


dL 

da 


- L . 




















An useful feature of the Hamiltonian H is the following well-known property: 

Proposition 2.5. Assume that the Lagrangian L in (|2.9I) does not depend on r (i.e., L = L(a, a, a)). 
Then H is constant along any solution a of (12.101) . 

The verification of the previous proposition amounts to a straightforward, direct computation 
of H along a, so we omit the details. In this case, we say that H is a prime integral of the 
biharmonicity equation. 


3. Biharmonic, rotationally symmetric, conformal maps between models of 

CONSTANT SECTIONAL CURVATURE 


The difficulty of the general problem (fourth order equations) suggests to restrict investigation to 
speficic, geometrically significant, families of maps. In particular, in this section we look for rota¬ 
tionally symmetric, proper biharmonic conformal diffeomorphisms between m-dimensional models 
of constant sectional curvature. We shall be able to obtain a complete description of such maps. 

First, let us point out when a rotationally symmetric map as in (12.41) is conformal. Comparing 
dilations of vectors which are respectively orthogonal and tangent to the radial direction, we easily 
find that a map of the type (|2.4I) is conformal iff 


d 2 


h 2 (a ) 

P{r) 


or, taking into account the boundary condition (12.51) . 


(3.1) 


h(a 

/ 0 ) 


Using (13.11) into the biharmonicity equation (12.111) . we find that the condition for biharmonicity of 
a conformal map of the type ([2.41) is the following: 


(m - T)f m 5 h(a) (/ 2 / (3) + h'(a) (4//" + (m - 5 )h{a)h" {a)) 

+ (3m — 14 )f' 2 h'{a) — 2 (m — 4)/ ,3 + 

/'((m — 7 )ff" — 2 (m — 4 )h(a)h"(a) — 2 (m — 4)h'(a) 2 ) — 
h(a) 2 P 3 \a) + (m — 2 )h , (a) 3 ' S j = 0 . 

In particular, if m = 4, equation (13.21) becomes: 

ff 2 f{ 3) + h\a){Aff" - h{a)h"(a)) 

(3.3) J V 

— 2 f 2 h\a) — 3 ff'f" — h(a) 2 h ly3 \a) + 2 h'(a) 3 ^j = 0 . 

Now, we are in the position to provide a complete description of solutions in the case of maps 
between 4-dimensional models of constant sectional curvature. More precisely, we have: 


Proposition 3.1. Let us consider rotationally symmetric maps as in (I2.4|) - (12.5I) . Let us assume 
that m = 4 and that both models have constant sectional curvature (see Remarks \2.1\ and IS.Sf) . 
Then the biharmonic conformal diffeomorphisms of type m-m can be enumerated as follows 
(c, d denote real positive constants): 

Case 1 - f(r) = r . 

A - h(a ) = a, a(r) = cr (harmonic diffeomorphisms from M 4 to itself); 

5 










B - h(a) = (1 /d) sin(cZa) and 


(3.4) 


(3.5) 


Case 2 


Case 3 


a (r) = — arctan(c 2 r), 

(proper biharmonic diffeomorphisms from M 4 to S 4 (d 2 ) \ {south pole}/; 
C - h(a) = (1/d) sinh(da) and 


2 1 

a(r) = — tanh _ 1 (c 2 r) (0 < r < -^) 

d c z 

(proper biharmonic diffeomorphisms from i? 4 (l/c 2 ) (Te., f/ie open ball of radius ( 1 /c 2 ) 
in M 4 / to H 4 (—d 2 )). 
f(r) = ( 1 /c) sin(cr). 

A - if h(a) = a, then there is no solution; 

B - if h(a) = (1/d) sin(da), then 


a(r) = (c/d)r (0 < r < —) 

gives rise to a harmonic conformal diffeomorphism, but in this case we do not have 
proper biharmonic examples; 

C - if h(a) = (1/d) sinh(da), then there is no solution. 
f(r) = ( 1 /c) sinh(cr). 

A - if h(a) = a, then there is no solution; 

B - if h(a) = (1/d) sin(da) ; then there is no solution; 

C - if h(a) = (1/d) sinh(da), then 


a(r) = (c/d)r (r > 0 ), 


produces a harmonic conformal diffeomorphism, while there is no proper biharmonic 
example. 


Proof. The proof amounts to a case by case analytical inspection of (13.31) . Indeed, 

Case 1 - f(r) = r. In this case, equation (13.31) becomes (here h = h(a)): 

2 h' 3 — h' (2 + h h") - h 2 h'" = 0 

which is satisfied for 

A - h(a) = a (harmonic maps); 

B - h(a) = (1/d) sin(da): in this case we have proper biharmonic diffeomorphisms whose explicit 
expression (13.41) for a(r) can be derived directly by the condition of conformality (I3.1IV ): 

C - h(a) = (1/d) sinh(rj a), which again produces proper biharmonic diffeomorphisms (also in this 
case, the explicit expression (13.51) for a(r) can be derived directly by the condition of conformality 
(13.11) 1. It is important to point out that, in this case, solutions are defined on a finite interval only. 
Case 2 - f(r) = (1/c) sin(cr). The biharmonicity condition (13.31) becomes 


(3.6) 2 h'(a) (cos(2cr) — h(a)h"(a) — 3) + 2sin(cr) sin(2cr) — 2 h(a) 2 h^(a) + Ah'(a) 3 = 0 . 

Therefore: 

A - if h(a) = a, the left-hand side of (13.61) becomes 



which cannot vanish identically, so that there is no solution in this case. 
B - if h(a) = (1/d) sin(da) (13.61) becomes 

4 sin 2 (c r) (cos (cr) — cos (da)) = 0 , 
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so that in this case we have harmonic diffeomorphisms, but not proper biharmonic examples. 

C - if h{a) = (1 /d) sinh(d a) the left-hand side of ()3.6I) becomes 

4sin 2 (cr)(cos(cr) — cosh(da)) . 

Now, since coshx > 1 for all x € R, it is obvious that the latter cannot vanish identically, thus 
there is no solution in this case too. 

Case 3 - f(r) = (1/c) sinh(cr). The biharmonicity condition becomes 

— 2 ^sinh(cr) sinh( 2 cr) + h^ 3 \a)h(a) 2 ^j 

(3.7) +2 h!(a) (cosh(2cr) — h(a)h"(a) — 3) + 4/t'(a ) 3 = 0 . 

Now, 

A - if h(a) = a, the left-hand side of (13.71) becomes 

— 8 sinh 2 ^ sinh 2 (cr) , 

from which we deduce nonexistence in this case. 

B - if h(a) = (1/d) sin(da), (13.71) becomes 

4sinh 2 (cr)(cos(da) — cosh(cr)) = 0 . 

As in Case 2C above, we easily conclude that there is no nontrivial solution. 

C - if h(a ) = (1/d) sinh(da), (|3.7I) becomes 

—4sinh 2 (cr)(cosh(cr) — cosh(da)) = 0 , 

which simply yields harmonic diffeomorphisms, but not proper biharmonic examples. □ 

Remark 3.2. We point out that (13.41) and (13.51) correspond to a dilation of R 4 composed with the 
inverse of the stereographic projection of 5 4 (d 2 ) and i7 4 (—d 2 ) respectively. 

Remark 3.3. A direct, case by case inspection shows that, if m > 2 and ra/4, then there exists 
no rotationally symmetric, conformal, proper biharmonic map between m-dimensional models of 
constant sectional curvature. By way of example, if f(r) = r and h(a ) = sin a, then condition (|3.2|) 
becomes: 

(3.8) 4 (m — 2) (m — 4) r m ~ 5 sin(2 a) sin 4 (a/2) = 0 , 

which is not possible when m ^ 2, 4. The other cases are similar, so we omit further details. 

Remark 3.4. It is interesting to compare the proper biharmonic diffeomorphisms from R 4 to 
S ' 4 \ {south pole} with the qualitative behaviour of rotationally symmetric harmonic maps (see 
El)- In particular, Jager and Kaul proved that the image of the functions a(r) associated to 
these harmonic maps cover a range [ 0 , R 4 ], with (vr/ 2 ) < < n, and a{r) oscillates around (vr/ 2 ) 

as r —> +00 (further details concerning the numerical value of R 4 can be found in [T3]). Jager 
and Kaul applied their results to draw some interesting conclusions concerning the existence of 
rotationally symmetric solutions to the Dirichlet problem for maps from the Euclidean unit m-ball 
B m to S m . In particular, if m = 4, they proved that the Dirichlet problem with boundary data 

(3.9) (0, 1) -»• (0, a( 1) = R*) 
admits a rotationally symmetric solution if and only if: 

(3.10) 0 < R* < R a . 

By contrast, the existence of the proper biharmonic conformal diffeomorphisms of Proposition 13.11 
implies that the boundary value problem (13.91) admits proper biharmonic solutions for all 

0 < R* < 7T . 
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Remark 3.5. The explicit solutions of Cases IB and 1C above were also studied in [3] and [ 161 . 
In particular, it was shown in m that they do not provide examples of biharmonic morphisms. 


Remark 3.6. If, in the context of Proposition 13.11 we admit solutions which are not continuous 
at the pole (i.e. which do not satisfy the boundary condition (12.51) 1. then we have some further 
examples of proper biharmonic conformal diffeomorphisms associated with the choice 


(3.11) 


a = — 


h{a) 

f ( r ) 


in place of (|3.1D . More precisely, by using similar arguments, it is not difficult to check that, in 
Case 1 ( f(r ) = r), we have: 

A - if h(a ) = a, then a(r) = c/r provides the proper biharmonic, conformal inversion p : 
R 4 \ {0} —>■ R 4 \ {0} given by p(x ) = cx/\x\ 2 . 

B - if h(a ) = (1/d) sin(da), then a(r ) = (2/d) arctan(c 2 /r) gives rise to a family of proper 
biharmonic conformal diffeomorphisms p : M 4 \{0} —>• 5' 4 (d 2 )\{two poles}, obtained by composing 
the inversion (we refer to A above) with the inverse of stereographic projection. Similarly, 

C if h(a ) = (1/d) sinh(da), then a(r ) = (2/d)tanh _ 1 (c 2 /r) gives rise to proper biharmonic 
conformal diffeomorphisms p : R 4 \ {R 4 (c 2 )} —> \ {pole}. 

In the remaining cases (/(?’) = (1/c) sin(cr) or /(r) = (1/c) sinh(cr)), an analysis similar to that 
in the proof of Proposition 13.11 shows that there are no examples of proper biharmonic conformal 
diffeomorphisms which satisfy (13.111) . 


4. Further developments: conformal, proper biharmonic maps from R m to an 

771-DIMENSIONAL model 


A natural, general development of the work which we carried out in Section [3] is to study when 
a conformal, rotationally symmetric map as in (12.4|) is proper biharmonic. 

Since the general case appears to be difficult, we start with the case that the domain model is 
the Euclidean space. In particular, using the Hamiltonian function H defined in (12.121) . we shall 
obtain a complete answer in the case that the domain model is R 4 . 

In order to simplify the analysis, it is convenient to perform the following change of variable: 

(4.1) r = e t , f £l, /3(t) = a(e t ) , 

so that the boundary condition (|2.5D becomes 

(4.2) lim f3(t) = 0 . 

t—y— oo 


Then the reduced bienergy functional (12. 8 |) takes the following form: 


(4.3) 


E 2(<Pf}) = \ f 

A JR 


(3 + (m — 2 ) p — [m — 1) h{(3) h! (/3) 


1 2 


e (m-A)t dt _ 


Also, we observe that, in terms of (3 , the conformality condition (13.11) becomes 


(4.4) p = h((3) . 

Here, and also in Section [5] below, we shall refer to the same rotationally symmetric map by writing 
ip a or pp, the choice for a. or /3 to describe the map being the one which is more suitable to simplify 
a specific discussion or calculation. Next, by using (12.111) . we compute explicitly the condition of 
biharmonicity. Together with conformality (14.41) and m = 4, that leads us again to the condition 
of Case 1 of the previous section, i.e., 

(4.5) 2 ti 3 - ti (2 + h h") - h 2 b!" = 0 . 


Now, we can state our main result: 












Theorem 4.1. Suppose that, in the notation of m, 


(4.6) 


: M 4 M' h \o') 

(9, r) ( 9 , a(r )) . 


is a rotationally symmetric, proper biharmonic conformal diffeomorphism. Then ip a is either as in 
Case IB or 1C of Provosition HOI 

Proof. Using the conformality hypothesis (14.41) and m = 4, we compute explicitly the Hamiltonian 
H (defined in (12.121) 1 associated to (14.31) : we obtain (up to a constant) 

(4.7) H = h 2 (1 — h' 2 + h h") . 

Since the Lagrangian L in (14.31) . when m = 4, does not depend on t, according to Proposition 12.51 
we conclude that H is constant along any solution. Now, since (14.21) holds and h( 0) = 0, h'{ 0) = 1, 
we deduce that the only possibility is that H = 0 along a solution. Therefore, in the presence of a 
nontrivial solution (i.e., /3 ^ 0) we deduce that necessarily 

(4.8) 1 - h' 2 + h h" = 0 . 

Next, derivation of (14.81) along a nonconstant solution leads us to conclude that 

(4.9) hh'"-h'h" = 0. 

Now, using the Jacobi equation (|2.3I) into (|4.9j) . we deduce that, along a nonconstant solution, 

(4.10) h'" = -Kh'. 

On the other hand, taking derivatives on both sides of the Jacobi equation (12.31) . we obtain 

(4.11) h!" = -Kh! -K'h . 

Finally, comparing (14.101) and ()4.11|) . we conclude that K' = 0, from which it follows that the 
solution is either as in IB of Proposition 13.11 (constant positive curvature) or as in 1C (constant 
negative curvature). □ 

Remark 4.2. In the final part of the proof, we have implicitly used the property that a model 
with constant radial curvature is a space form (see ED- 

4.1. Maps from a cylinder to a model. We show that the use of the Hamiltonian H can be 
useful in other related contexts. We consider rotationally symmetric maps as in (|2.6I) and assume 
that f(r) = 1, r £ M . The bienergy functional has (up to a constant) the following form: 

p+OO 

(4.12) 


1 /*—1~ OO 

E 2 ((p a ) = - / [a - X h(a)h'(a)] 2 dr. 

^ J — OO 


In this case the Hamiltonian H, computed by means of (12.121) . is given by: 


1 


H = — af a + - r a (2d - r Q 


(4.13) 

where r Q denotes the following expression: 

(4.14) r a = [a — A h(a)h’(a)] . 

Since L = L{a, a), according to Proposition 12.51 the Hamiltonian H is constant along solutions. By 
way of example, let us consider, in (14.121 ) . the special case of maps into a sphere (i.e., h(a) = sin ct). 
In [18] we pointed out that there are constant, proper solutions 

7T , 37T 


a 


; i n 1 j 


a = 


(4.15) 
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Here, by using the Hamiltonian H, we can prove that there are no other proper solutions (p a such 
that | T(ip a )\ = constant. Indeed, suppose that a gives rise to a nonconstant proper solution such 
that 

(4.16) |r(</> a )| 2 = [a — A sin(a) cos(a )] 2 = C 2 {C > 0) . 

Under these assumptions, the Hamiltonian H in (|4.13l) becomes 

(4.17) H = ±j(a-(±)C) , 

from which we easily deduce that a is a constant. It follows that the quantity sin(a) cos(a) must 
also be a constant, a fact from which our assertion follows immediately. By the same argument 
(simply use either h(a) = a or h(a ) = sinh(a) in (14.121) 1. we can also deduce that there exists no 
proper biharmonic critical point of (14.121) such that |T(y? a )| = constant, neither in M m nor in H m . 


5. Equivariant stability 


In this section we carry out the study of the second variation at the proper biharmonic con¬ 
formal diffeomorphisms which we described in Proposition 13.11 More precisely, we shall study the 
equivariant stability of the solutions given in (13.41) and (|3.5I) . Since we are only interested in the 
case that the domain model is M 4 (or the ball B 4 C M 4 ), it is convenient to make again the change 
of variable 61]). For the sake of convenience, we rewrite as follows the relevant bienergy (14.31) in 
terms of /3: 


(5.1) 


E 2(vp) = 


U 


P + 2/3-3 q((3) 


1 2 


dt 


where the integral is over the domain of (3 and, to simplify the calculations of this section, we have 
set 


(5.2) q(P) = h(P)h'(P) . 

Note that we shall only be concerned with the cases h((3) = (1/d) sin (d/3) (Case IB of Proposition 
EE} and h(/3) = (1/d) sinh(d/3) (Case 1C of Proposition 13.11) . For future use, we observe that, in 
terms of /?, the explicit expressions of the proper solutions (13.41) and (I3.5j) are 


(5.3) 
and 

(5.4) 


/3(t ) = - arctan(c 2 e 4 ) , t G 


(c> 0 ) 


P(t) = - tanh i (c 2 e t ) , t < In ( — ) (c > 0 ) 


respectively. We shall say that a critical point ipp of the type (12.41) is equivariant stable if 

d 2 E 2 ((p(i3 +S v)) 


(5.5) 


V 2 £ 2 (^)(V,V) = 


ds 2 


> 


0 


s=0 


for all smooth vector fields V = V(t) d/d/3, such that V(t) is a real valued, compactly supported 
function defined over the domain of /3(t) (V(t) ^ 0). 

In the following lemma we compute explicitly the equivariant second variation in (15.5} . 

Lemma 5.1. At the conformal, proper critical points of Proposition (Cases IB and 1C, i.e., 
63} and (15.41) respectively) we have: 

(5.6) V 2 E 2 (vp) (V, V) = J | [V + 2 V - 3 q'((3) v] 2 + 6 </'(/?) h(P) (//(/?) - 1 ) U 2 1 dt , 
where the integral is over the domain of f3 (t) . 
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Proof. First, we show that, for a generic critical point, 


(5.7) V 2 E 2 (vp) (V, V) = J { [v + 2 V - 3 q'(0) y] 2 - 3 q" (0) V 2 [/3 + 2 0 - 3 q({3) ] } dt 


To this purpose, the first step is to compute 
( 5g ) d_ E 2{P(p+sV)) _ 


k 

/ 


ds 

d_ 

ds 


0 + sV + 2(0 + sV) — 3 q(0 + si/) 


1 2 


dt 


0 + sV + 2(0 + sV)-3q(0 + sV) ■ V + 2V - 3 q'(/3 + sV) V 


dt . 


Next, using (15.81) : 

(5 9 ) f E ^VU3 +sV)) 


ds 2 
d 


ds - 


s =o 


+ sV + 2(/3 + sV) -3q(P + sV)) [V + 2V-3q'{(3 + sV)V 


8 = 0 


dt 


1 2 


V + 2V-3 q\/3)V -3 q"{0)V 2 0+2/3- 3q(0) 


dt . 


as required to prove (15.71) . Finally, we use the fact that our critical points are conformal and so 
they verify (14.41) . Taking derivatives, from (14.41) we deduce 

(5.10) (3 = h'(f3)/3 = h\P)h(f3) = q(/3) . 

Now, using (14.41) and (|5.10j) into (15.71) we obtain (15.61 ). as required to end the proof of the lemma. □ 


Remark 5.2. Integrating by parts (|5.7j) and using the fact that V(t) is a compactly supported 
function we obtain 

V 2 S 2 (^)(V,V) = J ( I(V),V)dt , 

where 

(5.11) /(V) = {V- [4 + 6q\(3)]V+ [9q ,2 (/3) - 3q" ((3) (f3 - 3q(f3))]v}. 

Direct computation shows that the operator /(V) defined in (|5.11|) can be derived by using the 
second variation formula obtained by Jiang in DUE] (see also {T9]). 


We are now in the right position to state our results in this context. 


Theorem 5.3. Let 

W : B\ 1/c 2 ) ->■ H\-d 2 ) 

be the rotationally symmetric, proper biharmonic conformal diffeomorphism defined by means of 
the function 0 in {53D ■ Then tpg is equivariant stable. 


Proof. By using h(0) = (1/d) sirih(d 0) into (|5.2I) and (15.61) . we hnd that, in this case, the second 
variation takes the following form: 


(5.12) X7 2 E 2 (vp)(V,V) = 


V + 2 V — 3cosh(2 d0) V +24 cosh(d 0) skilled/3) (cosh(d/3) — 1) V 2 > dt 


from which the conclusion of the proof is immediate. 

li 


□ 




























In a context of Dirichlet’s problem, we find that our conformal, proper biharmonic diffeomor- 
phisms into S ' 4 are equivariant stable if we restrict them in such a way that only the closed, upper 
hemisphere of S ' 4 is covered. More precisely, we have 


Theorem 5.4. Let us consider proper biharmonic solutions of the boundary value problem 


(5.13) 


pp : £ 4 ( 1/c 2 ) -> S 4 (d 2 ) 
fi(ln(l/c 2 )) = (n/(2d)) . 


The conformal, rotationally symmetric, proper biharmonic solution pp, defined by means of the 
function (3 in <Q with t < In ( 1 /c 2 ), is equivariant stable with respect to variations which preserve 
the boundary data. 


Proof. Formally, we have to consider the equivariant second variation (15.61) with respect to com¬ 
pactly supported vector fields V = V(t) d/d(3 which, in order to preserve the boundary data, satisfy 
the additional condition 

V (ln(l/c 2 )) = V (ln(l/c 2 )) = 0 . 

By using h(/3) = (1/d) sin(d/3) into (15.21) and (15.61) . we find that, in this case, the second variation 
takes the following form: 


(5.14) V 2 E 2 (pp)(V,V) = 

An(l/c 2 ) 


1 2 


V + 2 V — 3cos(2 d/3) V +24 cos(d/3) sin 2 (d/3) (1 — cos (d/3 )) V 2 > dt . 


Now, since 0 < f3 < (vr/(2d)) when t < ln(l/c 2 ), the conclusion of the proof follows easily from 
(15.141) . It should be noted that in this case the map pp covers exactly half sphere minus the 
pole. □ 


Remark 5.5. The methods of this section can be used to check that the proper biharmonic 
conformal diffeomorphisms described in A and C of Remark 13.61 are equivariant stable. Similarly, 
the proper biharmonic examples described in B of Remark 13.61 restricted to R 4 \ {B 4 (c 2 )}, are 
equivariant stable with respect to variations which preserve the boundary data (i.e., /3(lnc 2 ) = 
n/(2d)). 


Remark 5.6. We take this opportunity to point out that, in a similar stability context, formula 
(3.9) in [183 should be replaced by 


v- 


‘Epa’)(V,V)= / 
Jo 


p2tt 

2 V 2 - 2 V 2 k 4 

lo 

■ 


dd , 


from which it follows easily that the critical points of Theorem 3.2 in 
3.3 of [T 8 ] should then be deleted. 


are unstable, and Remark 
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